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When a microorganism begins swimming from rest in a Newtonian fluid such as water, it rapidly 
attains its steady-state swimming speed since changes in the velocity field spread quickly when the 
Reynolds number is small. However, swimming microorganisms are commonly found or studied 
in complex fluids. Because these fluids have long relaxation times, the time to attain the steady- 
state swimming speed can also be long. In this article we study the swimming startup problem in 
the simplest liquid crystalline fluid: a two-dimensional hexatic liquid crystal film. We study the 
dependence of startup time on anchoring strength and Ericksen number, which is the ratio of viscous 
to elastic stresses. For strong anchoring, the fluid flow starts up immediately but the liquid crystal 
field and swimming velocity attain their sinusoidal steady-state values after a time proportional to 
the relaxation time of the liquid crystal. When the Ericksen number is high, the behavior is the 
same as in the strong anchoring case for any anchoring strength. We also find that the startup 
time increases with the ratio of the rotational viscosity to the shear viscosity, and then ultimately 
saturates once the rotational viscosity is much greater than the shear viscosity. 


I. INTRODUCTION 

The hydrodynamics of natural and artificial micro¬ 
scopic swimmers in Newtonian and complex fluids con¬ 
tinues to be an active area of research [ 1-5] . Many of 
these studies focus on steady-state swimming. However, 
natural swimmers start, stop, and change direction, and 
artificial swimmers must do the same to be of use. There¬ 
fore, it is of interest to study transient swimming prob¬ 
lems, such as the acceleration of a swimmer from rest. 
In a Newtonian fluid of viscosity p and density p , the 
relevant time scale t v (V for ‘viscous’) for the startup 
of flow is the time t v = p/(pq 2 ) for changes in velocity 
to spread over a distance of order 1/q. For water and a 
length scale 1/q « 1 /im, we have t v ~ 1 ps, which is much 
shorter than the characteristic beat or rotation frequency 
of swimming microorganisms, such as 0.001 s for Vibrio 
alginolyticus [6], 0.01s for Escherichia coli [7], 0.02 s for 
Chlamydomonas reinhardtii [8] and sea urchin sperm [9] , 
0.05 s for human sperm [10], and 0.5 s for Caenorhab- 
ditis elegans [1 1]. Calculations for how the swimming 
speed rapidly rises to its steady state value for an ide¬ 
alized swimmer in a Newtonian fluid were presented by 
Pak and Lauga [12], Complex fluids have additional time 
scales, which are much longer than t v , and can be com¬ 
parable to or even longer than the characteristic beat pe¬ 
riod. For example, polymer solutions in which swimming 
has been studied can have relaxation times of the order 
of seconds [11, 13]. Likewise, the airway mucus encoun¬ 
tered by beating cilia has relaxation times of the order 
of tens of seconds [14, 15]. The startup problem for a 
swimmer in a viscoelastic fluid has been examined by El¬ 
fring and Lauga [16]. Swimming bacteria have also been 
recently studied in liquid crystal solutions [17, 1 ]. For a 
liquid crystal with Franck elastic constant K ss 10 pN, 
and shear viscosity p ~ 10 Pa-s [17], the characteris¬ 


tic relaxation time t e (‘e’ for ‘elastic’) for distortions 
of the liquid crystal with length scale 1/q ~ 1 pm is 
t e = p/(Kq 2 ) ss 10 s. The long relaxation times of these 
complex fluids can lead to much longer-lived transient 
swimming flows than in the Newtonian case. 

In this article we study the startup problem for an 
idealized swimmer with small amplitude waves in a hex¬ 
atic liquid crystal Him. The hexatic phase is studied 
because its theory is mathematically simpler than that 
of the nematic phase, yet it retains some of the same 
distinctive features such as Franck elasticity, rotational 
viscosity, and anchoring effects. We show that the time 
required to attain the steady swimming speed is propor¬ 
tional to the relaxation time of the liquid crystal. Since 
this time is much longer than the viscous startup time 
t v , the swimmer first attains the Newtonian swimming 
speed in a time comparable to t v , and then reaches the 
final swimming speed over a longer time t e . We study 
how the evolution of the swimming velocity depends on 
the anchoring strength of the the liquid crystal at the sur¬ 
face of the swimmer, and also the Ericksen number Er, 
which is the ratio of viscous to elastic stresses. For strong 
anchoring, we find that the swimming speed reaches its 
ultimate value after a time proportional to the Ericksen 
number. At high Ericksen number, the behavior of the 
flow field, the liquid crystal configuration, and the swim¬ 
ming speed is independent of anchoring strength, and is 
given by the strong anchoring case. When the anchor¬ 
ing strength vanishes, we find that the swimming speed 
can oscillate as it approaches its final value, and that the 
swimmer can even reverse direction several times before 
reaching its steady speed. 

The outline of the paper is as follows. Section II de¬ 
scribes the governing equations and timescales. Then 
the case of strong anchoring is solved analytically in 
Sec. III. Section IV considers the case of arbitrary an- 
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Figure 1. (Color online) Sketch of the one-dimensional swim¬ 
mer and definition of the angle field 6. The liquid crystal 
particles are not drawn to scale. The outward pointing unit 
normal vector on the swimmer is N, and the tangent vector 
to the swimmer is l. 

choring strength in the limits of high and low Ericksen 
number, paying special attention to the striking case of 
zero anchoring strength. The conclusion is Sec. V. 


II. SETUP AND GOVERNING EQUATIONS 

A. The swimmer 

The swimmer is a line with a prescribed traveling wave. 
In the frame of the swimmer, the material points on the 
swimmer have coordinates 

x m = x ( 1 ) 

y m = bsm[qx-u](t)t], ( 2 ) 

where b is the amplitude of the wave, q is the wavenum¬ 
ber, and w(f) is the time-dependent frequency (Fig. 1). 
Following Taylor [19], we assume the amplitude is small 
compared to the wavelength, b 2 q 2 <C 1. Initially, the fre¬ 
quency vanishes, w( 0 ) = 0 , and the fluid is motionless. 
The frequency rises to its final value io Q0 with a character¬ 
istic time scale t s (‘s’ for swimmer). This characteristic 
time is long compared to the viscous startup time t v , but 
seems to be comparable to or longer than the beat pe¬ 
riod 27r/w 00 , depending on the species. For example, the 
switching time from backward to forward swimming, or 
vice versa, in V. alginolyticus is of the order of tenths of 
seconds [ 20 ], about ten times longer than the beat period. 
When sea urchin spermatozoa are rendered immotile by 
lowering the pH, the waveform takes a few periods to 
restart once the pH is raised again [21]. 

B. Hexatic Liquid Crystals 

The swimmer lies in a two-dimensional hexatic liquid 
crystal. Recall that a hexatic liquid crystal has six-fold 
bond-orientational order, described by an angle field 9, 
where 9 and 9 + 2tt/6 denote the same physical configu¬ 
ration [22] (Fig. 1). We use the notation and conventions 
of reference [23], and refer the reader there for more de¬ 
tails on the derivation of the governing equations. We 


consider incompressible flow, V ■ v = 0, where v is the 
velocity. The dynamical equations for velocity and angle 
field are 

-V P + pV 2 v - K{V9)S7 2 9 + ^ V x (zV 2 6») = 0(3) 

dtO + v • V0 - ■ V x v = —V 2 9,( 4) 

2 7 

where p is shear viscosity, K is the Franck elastic con¬ 
stant, z is the unit vector perpendicular to the plane 
of the liquid crystal film, and 7 is the rotational viscos¬ 
ity [24] . Note that the pressure p has been defined so that 
it vanishes in equilibrium [23]; in other words, p is the ac¬ 
tual pressure minus the equilibrium pressure arising from 
Franck elasticity [22], Note also that we have assumed 
the Reynolds number to be vanishingly small. This ap¬ 
proximation is valid because the liquid crystal solutions 
which motivate our study are 10000 times as viscous as 
water [17], and also because we do not attempt to resolve 
the dynamics on the scale t v . 

At the swimmer, we assume no-slip boundary condi¬ 
tions on the velocity field: 

v(z m ,t/ m ,t) = y m y. (5) 

Far from the swimmer, the velocity field has an unknown 
uniform value v(j/ —> 00 ) = f/x that we must solve for. 
The flow U in the rest frame is the swimming velocity 
in the lab frame, which has opposite sign to U. We 
also impose anchoring conditions on the liquid crystal 
near the swimmer surface. Since we assume the swimmer 
has a small slope, we may assume the angle <j> between 
the swimmer tangent vector £ and the ir-axis is small. 
With accuracy of order b 2 q 2 , the anchoring condition be¬ 
comes [2 ] 

A'N • V9 + W{9 - <t>) = 0, ( 6 ) 

where N is the outward-pointing normal vector and W 
is an energy per unit length giving the strength of the 
anchoring potential. 

As mentioned earlier, the initial condition for the fluid 
is v = 0. The initial condition for the angle field is the 
equilibrium configuration, which satisfies V 2 0 = 0 and 
the anchoring condition Eq. ( 6 ). 


C. Time scales, nondimensionalization, and the 
impulsive startup problem 

There are four time scales in our problem, with a strong 
separation of the viscous time scale from the others: 

t v <C w ” 1 <t s < t e . (7) 

To make the governing equations dimensionless, we use 
l/cuoo as the units of time and 1/q as the units of length. 
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Then, with p measured in units of yuioo, we have [ 

- Vp + V 2 v - V6(y 2 0) + -^VV 2 « x z = 0 ( 8 ) 

rL/r 2hji 

d t 8 + v • V0 — -z • V x v — - — V 2 0 = 0,(9) 
2 7 Er 

where 


Er = 


Ahf ' 


( 10 ) 


It is convenient to use the stream function ip, which auto¬ 
matically enforces the constraint of incompressibility by 
the definition v = V x (ipz). In terms of ip, the governing 
equations are 

VV + -^V 4 0 + 2-z ■ VO x W 2 0 = 0 (11) 

2hji hr 

dtO + lV 2 ip + Z • VO x Vip - ^ 2 -' V 2 0 = 0. (12) 
2 7 Er 

In dimensionless form, the boundary condition for the 
angle field 0(x m ,y m ) at the swimmer is 


N • VO + w(8 — </>)= 0, (13) 

where w = Wj ( Kq ). The dimensionless no-slip boundary 
condition at is 

. . d\oj(t)t] . , . , „ . 

v(x m ,y m ) = -e———cos [x — uj(t)t]y, (14) 

where we have introduced e = bq. Note that we are 
also measuring frequency in units of Woo To simplify the 
analysis, we will consider the impulsive startup problem 
and suppose that uJooT <C 1 , even though r is typically a 
longer time scale than I/wqq. Therefore, we suppose the 
dimensionless frequency impulsively jumps from zero to 
unity at t = 0. Denoting by H(t) the unit Heaviside step 
function, we have 

v(x m ,y m ) = cos [x - H(t)t] y. (15) 

Our task is to solve for the swimming speed —U as a 
function of time. We expand in powers of e, 

V = v ( 1 ) +v (2) + --- (16) 

0 = 6> (1) + 6> (2) + • • • (17) 

where the superscript denotes the power of e. 


III. STRONG ANCHORING 

First we consider the case of strong anchoring, w —> oo, 
since in this limit we can get explicit expressions for 0 
and U as functions of time for impulsive startup of the 
beating. 

In terms of the stream function, the first-order parts 


of equations ( 11 - 12 ) are 


V 4 ?/> (1) + JrV 4 fl (1) = 0 

2 br 

(18) 

d t 0 {1) + - v 2 ^, (1) - -^-V 2 0 (1) = 0 . 

2 7 E 1 ' 

(19) 

The initial conditions are 


0 (1 ) (x, y, t = 0 ) = ee~ v cos x 

( 20 ) 

ip (1) {x,y,t = 0 ) = 0 , 

( 21 ) 

and the boundary conditions are 


V^ (1) \ y=0 = -e^sin [x-tH(t)]x 

( 22 ) 

0^\ y= 0 = ecos(x — t). 

(23) 


Note that the angle field has some initial distortion due to 
the strong-anchoring condition. Furthermore, the initial 
condition for the angle field has exactly the same spatial 
form as the sinusoidal-steady state solution for 0 ^ that 
was found in [2. ]. Thus, there is no need for distortions 
in the angle-field to spread once the swimmer waveform 
starts up; the angle held simply starts to oscillate in time 
with the same frequency as the swimmer waveform. The 
fluid velocity, however, changes discontinuously, since we 
do not resolve dynamics on the small timescale t v (stud¬ 
ied in [12]). For t > 0, we find 

v^ = —eye~ v sin(x — t) (24) 

= —e(l + y)e~ v cos(x — t) (25) 

O^ = ee~ v cos(x — t). (26) 

These results can also be found by directly solving the 
equations (18-19) using Laplace transforms. 


To see how the final sinusoidal steady-state solution 
emerges, we must turn to the second-order equations: 


- Vp (2) + vV 2) 


1 

Er 


-V x (zV 2 6 » (2) ) - V 6» (1 W 1} 


(27) 


and 


— ^V 2 6 » (2) =<% 6 > ( 2 ) +vW-V 6 >W — -z-Vxv< 2 >. (28) 
Er 7 2 


Averaging over a spatial wavelength, these equations be¬ 
come 

d 2 y (i 4 2) > + 4r^ (0(2)) = f (29) 

= ^ 2 < 0 (2) > - \{d y vM) - g , ( 30 ) 

where / = (d x 0^ V 2 0^)/Er and g = (vW • VfH 1 )). For 
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Figure 2. (Color online) Time evolution of the average second- 
order angle field for Er = 1 and /r = 7, or D = 5/4. 

Starting from the top of the figure, the first eight curves (blue) 
correspond to u>oot = 10 —2 , 10 1 , 10 5 , and the last curve 

(black) corresponds to the steady-state value, or w ™t —¥ 00 . 


strong anchoring, 

/ = 0 (31) 

g = (e 2 /2)(l + 2y)exp(-2y)H(t). (32) 

The initial conditions are (0( 2 ))|t=o = 0 and (v^ 2 ^)| t= o = 
0, and the boundary conditions are 

(4 2) )| y= o = -(VmdyV^ly^o = e 2 /2 (33) 

(9^)\ y=o = 0. (34) 

Note that to arrive at (33-34), we expanded the bound¬ 
ary condition at y = y m and used the fact that the 
strong-anchoring angle field is out of phase with the wave 

y = y m - 

Using C to denote the Laplace transform, 

/»00 

£{/(*)} = / exp(—st)f(t)dt, (35) 

Jo 


let V(s) = £(v^), and 0(s) = C(9(‘ 2 '>). Then 


S0 


V" 
Er 7 


2Er 


0 '" = 0 


0 " 


y = -9, 


(36) 

(37) 


where g = g/s and the primes denote differentiation with 
respect to y. Note that (36) and the fact that 9^ and 
Vx ' are finite at y — >• 00 implies that V' + 0"/2/Er is a 
constant. Thus, (37) becomes diffusion with a source, 


s0 - DQ" = -g, 


where the diffusion constant is 




(38) 


(39) 


§ _ e 2 a - 12£> 

2s (s - 4 D) 2 

e 2 12 D — s + (8 D — 2 s)y _ 2y 


2s 


(s - ADf 


(40) 


Note that we recover the proper long-time limit of the 
second-order angle field in steady-state swimming [23]: 


iST 6 = ‘ 2 2(4^7) [O + ^e-^-a] ■ (41) 


The angle field 0 has an explicit inverse Laplace trans¬ 
form: 


(0 (2) ) = {( 6 + 4y ) e 2y - 6erfc 2//( 2 V A Dt) 

+ (8 Dt -3-2 y) 

x [e 4£>t_2y erfc (z-) — e 2y erfc (2+)] } , (42) 

where erfc is the complementary error function [25], and 


4 Dt ± y 
2 y/Dt ‘ 


(43) 


Note that since D oc 1/Er, we have (9^) oc Er. Figure 2 
shows how ( 6 ( 2 )) evolves to its final steady-state value. 
Note that for any finite t. the second-order angle field 
vanishes at y —» 00 . However, for a given y ft 1, we 
can always wait long enough for the angle field to reach 
its steady-state value of — 3e 2r yEr/(8y + 7). Despite the 
complicated form of Eq. (42), close examination of the 
plots of ( 9 ( 2 -*) for various t (Fig. 2) reveals that the y- 
value at which |(0^)| attains half its maximum value 
increases like t 1 / 2 , as expected for diffusion. 


We can also solve for the velocity field in Laplace space: 


V = + v 2 e~ 2y + U, (44) 


where 

2 y/s 3 /D - \2sfsD 

V-\ = 6 - 

4Ers(s — AD) 2 
2 — Asy + 16D(1 + y) 
V2 = 6 4Ers(s - 4L>) 2 ’ 


(45) 

(46) 


and the term U remaining in the limit y —> 00 is the 
swimming speed in the frequency domain, 


J_ _ 4 + y/s/D 

2s sEr(2/D + y/s ) 2 


(47) 


Note that the long-time limit of the swimming speed 
agrees with the steady-state result [23] : 


lim sU = Uao 

s —^0 


e 2 Ay — 7 
2 Ay + 7 


The solution is 


(48) 
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Figure 3. (Color online) Dimensionless swimming speed vs 
dimensionless time on a log-log scale for fj, = 7 and Er ranging 
from (from left to right) Er= 0.01 (blue), Er=0.1 (gold), Er=l 
(green), Er= 10 (red), and Er= 100 (purple). 


is a transient contribution that decays for large t as 


Ut r ---(54) 

16 (n /7 + l/4) 3/2 V nt 

The evolution of the swimming speed in time for vari¬ 
ous Er is shown in Fig. 3. Since the curves in Fig. 3 are 
equally spaced in log(Er) and are also equally spaced on 
the log-log plot, the time to reach the midway point be¬ 
tween uiooqb 2 /2 and Uoo increases linearly with Er. We 
can also see that the time to reach the steady state in¬ 
creases linearly with Er from Eq. (54): the dimensionless 
time t* it takes C4 r to decay to a small fraction of the 
steady-state speed is 


t* oc 


Er 


(p/t + V 4 ) 3 ' 


(55) 


Note that t* oc Er(y//i) 3 for 7 <C n, and t* oc Er for 
7 » + 


By inspection of 0, V, and U, we observe that the char¬ 
acteristic time scale for disturbances to diffuse over a 
distance L is L/y/D oc Ly/'Eii'y/lfj, + 7/4). 


We were unable to find an analytic form for the ve¬ 
locity field in the time domain. However, we can solve 
for the swimming speed in the time domain by inte¬ 
grating Eq. (36) in the time domain with respect to y. 
Let V = (vx 2 ^) and 0 = (O^). One integration yields 
V' + 0"/(2Er) = 0, since all y-derivatives of V and 0 
must vanish at y —> 00 . One more integration shows 
that V + 0'/(2Er) is constant; this constant must equal 
the flow speed at y —> 00, again because 0' vanishes at 
y —>• 00. Thus, 


U = V{y,t)+^ a -^(v,t). 


But V(0,t) = e 2 /2; therefore 

u = i + ^ 6 


y =0 


Substituting Eq. (42) into (50) yields 

U = Uoo + Ut r , 

where 

tj = I _ 1 = 1 7 

00 2 4DEr 2 4y + 7 


(49) 

(50) 

(51) 

(52) 


is the steady-state speed at t —> 00 [23], and 


U tI = e 


2 


1 nr 1 

2Er V Dtt Er 


—— t] e 4Dt erfc f 2 yfDt) 
AD J \ J 

( 53 ) 


IV. ARBITRARY ANCHORING STRENGTH 

We now turn to the swimming problem with finite but 
arbitrary anchoring strength w. Again we work in pow¬ 
ers of the dimensionless amplitude e. To first order we 
can find an analytic form for the Laplace transforms of 
the angle field and stream function. At second order we 
consider limits such as large Ericksen number, small Er- 
icksen number, and vanishing anchoring strength. 


A. First-order swimming problem 


To write the Laplace transform of the governing equa¬ 
tions (8-9) to first order in e, we need to find the ini¬ 
tial value of the angle field, since C{dtd^} = —9^(t = 
0) + The initial condition for the angle field is 

that it is the equilibrium field obeying the anchoring con¬ 
dition (6) for the initial swimmer shape, y m = —iexpia; 
in complex notation, which leads to 


0 = V 4 i^ (1) 


2Er 

s(j {1) = —V 2 ^ (1) - 
yEr 


y 4 0« 

l^ (1) 


The solutions to Eqns. (56-57) are 


we 

1 + w 


e ix ~ y 


(56) 

(57) 


0(D 


(c 0 + Cl! ,)e-^-^^ 


— Cl + 


■ w 


c 2 e 


ky 


(58) 

(59) 


where 


k 



4syEr 
7 + 4/it’ 


(60) 
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and the s-dependent coefficients Co, c±, and C2 are deter¬ 
mined by the no-slip and anchoring boundary conditions 
at the surface of the swimmer. In Laplace space, the no¬ 
slip boundary condition v(a; m , y m ) = (x m , dty m ) becomes 


(0„lP>,-0 x ip>)| 



i + s 


(61) 


to first order in dimensionless form. Likewise, to first 
order, the anchoring boundary condition (13) becomes 


dy~e {1) 


2/=0 


+ ID 


0(L 


— e- - 

y=0 1 + S 


= 0. 


(62) 


Although analytic expressions for Co, ci, and C2 and also 
i/K 1 ) and 9^ can be computed for arbitrary w, their 
lengthy form prevents us from displaying them here. In 
the next two subsections we consider the limits of large 
and small Ericksen number. Just as in the strong anchor¬ 
ing case and the sinusoidal steady-state case [23], there 
is no swimming speed to first order. 


B. Asymptotic solution at high Ericksen number 


Our analysis in the limit of Er 1 is similar to that in 
the sinusoidal steady-state problem [23]. Expanding the 
coefficients from the first-order problem to leading order 
in inverse powers of Ericksen number, we find 


Co = 


Cl = 


le 

i + s 
ie 

i + s 


+ 0| s 


C2 = 4= i,vl/4+ g^+oiJ- 


a/Et (i + s)\/s3 


Er 


(63) 

(64) 

(65) 


With an accuracy of 0(l/y/Et), the coefficients co, Ci, 
and C2 are independent of the anchoring strength, and 
take the same values as in the strong-anchoring problem; 
the rapidly varying component, proportional to expfcy, 
which does not arise in the strong-anchoring limit, has an 
amplitude of 0(l/vEr). Thus we conclude that the first- 
order flow field and angle field at large Ericksen number 
and arbitrary anchoring is well-approximated by the so¬ 
lutions to the strong- anchoring problem. 

The second-order flow and angle field at large Erick¬ 
sen number and arbitrary anchoring strength is also given 
by the strong-anchoring solutions. To see why, note that 
large Er, the second-order equations (29-30) are identi¬ 
cal to the strong-anchoring equations since / and g de¬ 
pend only on first-order quantities. For arbitrary w, the 
boundary conditions at second order are 


( 4 2) ) = ~(ym d v v^) ( 66 ) 


and 

{-d y 9^+w9 ^) y= 0 

= (~d x y m d x 9 w +y m dy9 (1) - wy m d y 9 {1) ) v=0 . (67) 

The terms on the right-hand sides of (66-67) arise from 
expanding the boundary about y = 0, and depend on 
the anchoring strength w implicitly through the first- 
order quantities. However, for large Er these quanti¬ 
ties are given by the strong-anchoring limit; in partic¬ 
ular, all the terms on the right-hand side of (67) vanish 
since y m and 9 ^ are out of phase for strong anchoring. 
Furthermore, the value of d y 9^ at y = 0 is O(Vfti). 
Since the angle field 9^ is 0(Er) [see Eq. (42)], the term 
(d y 9^) y= 0 is subleading in Er. Therefore, large Er be¬ 
havior of the arbitrary anchoring strength problem is the 
strong-anchoring problem. The steady-state angle field 
evolves slowly, with a characteristic time proportional to 
Er. 


C. Asymptotic solution for low Ericksen number 

In this section we show explicitly for arbitrary dimen¬ 
sionless anchoring strength w that the angle field attains 
its steady-state value quickly when the Ericksen number 
is low, with a startup time proportional to Er. Consider 
the governing equations (3) and (4) in the limit of Er -Cl. 
In this limit the viscous stresses are weak compared to 
elastic stresses. Note that sending Er to zero removes the 
highest (and only) time derivative in Eq. (4). Therefore, 
we have a singular perturbation theory problem, with an 
inner problem for early time and an outer problem for 
late time. The outer problem is the steady-state prob¬ 
lem studied in [23]. We introduce the short time scale 
r = fEr to resolve the fast relaxation of the angle field 
and flow to their steady-state values. Thus, we formally 
expand the fields at each order of e in powers of Er. For 
example, Er) = t ) + Er^^f, r) + .... Note 

that 


£> t 9 w = + d t 9 « + 8 T 9 + .... ( 68 ) 

Only the angle field enters the zeroth order equations: 

V 4 ^ = 0 (69) 

d T 9 { 0 1] = -V 2 9 ( k ) . (70) 

7 

Expanding the boundary condition in Er leads to 

w(^ 1) | !/= o - ecos(:c)) = d y 9^ ) \ y=0 . (71) 

To zeroth order in Er, the solution is 

= e _—— e ~ v cosx. (72) 

Note that since the wave of the swimmer is stationary 
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Figure 4. (Color online) Plots of the dimensionless swimming speed U/ {uJooqb 2 ) vs. r for Er <C 1 and log( 7 / y) = 0 (left panel), 
log( 7 /p) = 1 (middle panel), and log( 7 / y) = 3 (right panel). Colors correspond to w = 0.001 (black), w = 1 (purple), and 
w = 1000 (blue). Red dashed lines indicate the steady swimming speed found in [23]. 


on the fast time scale r, the angle field is simply the 
equilibrium angle field for a stationary ripple. 


Turning now to the next order in Er, we have 


To first order in Er, we have 

d 2 y (v$) + \dl{ef ] ) = 

dr{ef ] ) + \d y {v%) = ±d 2 y (e [ 2) > - (v^ ■ ve^).(79) 


3 ( 1 ) 


i 


V 4 ^ = 0 

ZnlX 1 ) — ^Y 7 


d T e\ L> + 

2 7 

with boundary conditions 

{d y ^\ - (0,-ecosa;) = 0 

w(9^ — ersina;) — d y 6 ^ 


= 0 


J y=0 


(73) 

(74) 

(75) 

(76) 


and initial conditions i/jq 1 ' 1 = 0 and 9^ = 0. The form 
of Eqs. (73-74) is similar to the form of (56-57), and the 
Laplace transforms of ^q 1 '* and 9 j 1 '* are readily found. 
However, we are not able to find an explicitly analytic ex¬ 
pression for the inverse Laplace transform of these quan¬ 
tities. The boundary conditions (75-76) and the govern¬ 
ing equations (73-74) imply oc sin a; and $J 1 ' ) oc sin a;. 
These phase relations apply because the swimmer wave 
is stationary on the small timescale fEr over which the 
fields develop from their initial values. 


The boundary conditions to second order in e and first 
order in Er are greatly simplified by the observation that 
9 1 and y m are in phase in x , which makes many terms 
vanish when we expand y rn to first order in Er, y m = 
e sin x — eErr cos x. Thus 

( v xo) \y=o=-e{d y v { JJ sin x) (80) 

[-(9^ 2) )+^(^ 2) )] y= o=° ( 81 ) 

Since is independent of time, it is straightforward to 
analytically determine the Laplace transform Eq. (79). 
Using the boundary condition (81) to determine the 
Laplace transform of ( 9\ ), we may use the same meth¬ 
ods as in Section III to solve for the swimming speed 
U(s). We could not find an exact inverse Laplace trans¬ 
form for this function except in the limit w —> oo, so 
we invert the expression numerically using a Fourier/de 
Hoog type method [26] on the interval t/r € [0.05,20] 
with 2,048 gridpoints. Figure 4 shows that the swim¬ 
ming speed relaxes to the steady state value with a char¬ 
acteristic time proportional to Er for various anchoring 
strengths. 


Now consider the problem to second order in e. To 
zeroth order in Er, the second-order problem is 

dy(do ) ) = 0 

dX 2) = ^$ 2) >, (77) 

(n\ 

suggesting that (9q ’) is a constant. The boundary con¬ 
dition at this order is 

[H^o 2) >+™<0o 2) >] y=o = O, (78) 

implying (9 { Q 2) ) = 0. 


D. Zero anchoring 

The final limit we consider is zero anchoring strength, 
vj = 0, for arbitrary Ericksen number. We proceed as 
above, by solving the governing equations to first order 
in e analytically using the Laplace transform, finding the 
angle field and swimming velocity in the frequency do¬ 
main to second order in e, and then using the Fourier/de 
Hoog method to invert the Laplace transform to find the 
swimming speed U(t). The results are shown in Fig. 5 for 
Er=l. When y/y is sufficiently small, less than around 
10 , the ultimate swimming direction is opposite the direc- 
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Figure 5. (Color online) Log-log plot of the dimensionless time 
t* to attain the steady-state swimming speed versus 7 / p. for 
Er = 1 . The time t* is defined to be the time in which the 
amplitude of the envelop of the oscillations (if there are any) 
is 5% of the ultimate swimming speed. The red dots are the 
values of t* for zero anchoring strength, w = 0. For com¬ 
parison, the solid blue line corresponds to the case of strong 
anchoring. To best compare the dependence of 7 /p, the time- 
origins of both graphs have been shifted by defining t* = 1 for 
7 — fi. The insets show the numerically computed swimming 
speeds versus dimensionless time for 7 / p = 1 , 10 , and 1000 . 


tion of the traveling waves of the stroke [23]. When 7/p is 
large enough, the swimmer ultimately swims in the same 
direction as the traveling waves. Figure 5 shows that in 
these cases the swimming speed oscillates about the final 
swimming speed as the steady state develops; when 7/p 
is large there can be several time intervals in which the 


swimmer changes direction before it settles down into its 
steady swimming speed. 

V. CONCLUSION 

In conclusion, we have studied the time-evolution of 
a swimmer in a hexatic liquid crystal him. The swim¬ 
mer in our approximation has a small amplitude sinu¬ 
soidal stroke that begins abruptly. Since the startup of 
viscous how is much faster than the evolution of the liq¬ 
uid crystal configuration, the swimmer immediately be¬ 
gins swimming with the Newtonian swimming speed. It 
reaches the ultimate steady-state swimming speed with 
a dimensionless characteristic time that is proportional 
to the Ericksen number, or, equivalently, a characteris¬ 
tic time proportional to p/(Kq 1 2 3 4 5 6 7 ). When the Ericksen 
number is large, the behavior is independent of that an¬ 
choring strength and given by the limit of inhnite anchor¬ 
ing. When the anchoring strength vanishes, the transient 
swimming speed can oscillate, and even change sign. Al¬ 
though the high symmetry of a hexatic liquid crystal 
lead to great simplifications in our analysis, we expect 
that many of the phenomena we found here will serve 
as a guide for the investigation of more the more real¬ 
istic but complicated case of a finite-size swimmer in a 
three-dimensional nematic liquid crystal. 
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